Classical diamagnetism, magnetic interaction energies, and repulsive forces in 

magnetized plasmas 
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The Bohr-van Leeuwen theorem is often summarized as saying that there is no classical magnetic 
susceptibility, in particular no diamagnetism. This is seriously misleading. The theorem assumes 
position dependent interactions but this is not required by classical physics. Since the work of 
Darwin in 1920 it has been known that the magnetism due to classical charged point particles 
can only be described by allowing velocity dependent interactions in the Lagrangian. Legendre 
transformation to an approximate Hamiltonian can give an estimate of the Darwin diamagnetism 
for a system of charged point particles. Comparison with experiment, however, requires knowledge 
of the number of classically behaving electrons in the sample. A new repulsive effective many-body 
force, which should be relevant in plasmas, is predicted by the Hamiltonian. 

PACS numbers: 75.20.-g, 52.27.Aj, 05.20.Jj 
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INTRODUCTION 



The Bohr-van Leeuwen [IH3 (BvL) theorem states 
that the magnetic susceptibility of a classical system of 
charged point particles interacting via some position de- 
pendent potential energy is zero. Further support for 
this theorem can be found in ■ Landau showed that 
quantum mechanics can explain diamagnetism in metals 
[l0(. Recent work, however, shows that perfect conduc- 
tors exhibit classical perfect diamagnetism 
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seem- 



ingly in blatant conflict with Bohr and van Leeuwen. The 



BvL theorem has also been questioned by Dubrovskii [13 1 
on the grounds that it neglects a relevant constant of the 
motion other than the energy. Usually the message of the 
BvL theorem is summarized as proving the nonexistence 
of classical diamagnetism Hij . We will show here 

that this is highly misleading by displaying an accurate 
energy expression for a system of charged point parti- 
cles in an external magnetic field. We then also discuss 
the Hamiltonian corresponding to this energy and draw 
some general conclusions about the behavior of magne- 
tized plasmas. 

It is obvious that one can treat magnetic susceptibility 
using classical models if one gives up the assumption of 
point particles. Using classical objects that are extended 
balls of charge one can find a classical explanation of dia- 
magnetism |16| , and classical models with dipoles can ex- 
plain paramagnetism. One can, however, reasonably ar- 
gue that these are not fundamental in the same way that 
point monopole particles are. Since the work of Charles 
Galton Darwin [17] in 1920 it has been known that the 
correct Lagrangian for a system of charged point parti- 
cles requires velocity dependent interactions. Once the 
BvL assumption of only position dependent interactions 
is relaxed one finds classical diamagnetism effortlessly, as 
we now proceed to show. 



ON THE MAGNETIC ENERGY OF A SYSTEM 
OF POINT CHARGES 

The electromagnetic energy of a system can be ex- 
pressed in many ways. Here we first assume that we 
are dealing with charged point particles in vacuum, i.e. 
no dipoles. It is then sufficient to consider the two fields 
E and B, or equivalently the potentials <fi, A. The var- 
ious expressions for the energy of a system, and their 
interrelations, that then can be written down have been 
reviewed by Franklin 18| . One well known expression for 
the energy is given by, 
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The electric energy is however normally taken into ac- 
count by finding the electrostatic potential, (f>(rj\rk) at 
particle j, due to the other particles k ^ j of the system. 
This gives, 



where self-interactions are assumed removed. Consider 
now the magnetic energy. A moving charged particle 
produces the magnetic field, 



B 3 (r) 



e 3 Vj x (r - r s ) 



c r — r.,- 



(3) 



to first order in v/c. The total (internal) field is then 
B- l {r)=Y,B 3 {r) 



(4) 



To estimate the energy we should then introduce this 
in the integral and integrate over all of space. To get 
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finite results one must again ignore self-interactions. This 
means that we put, 
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(5) 



for the internal magnetic interaction energy. The calcu- 
lation gives (Breitenberger |19j|). 

\- e e k [v r v k + (v r e k j){vj-e kj )} 



j<k 
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Here r kj = |r 3 - - r k \ and e kj = (r,- - r k )/r kj . This 
is the Darwin magnetic energy expression that follows 
without approximation from the Darwin Lagrangian. By 
introducing the internal vector potential, 



Ai(rj;r k ,v k ) = ^ 



e k v k + (v k ■ e kj )e kj 
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one finds, 



V — 

3 



^ 2c 



Vj ■ Ai(rj;r k ,v k ), 



(7) 



(8) 



as an alternative expression. 

If there is also an external magnetic field, B c (r), the 
magnetic energy will be, 



J (Bi + B c ) 2 dV = E ini + E ic + E c 



(9) 



This form of the magnetic energy makes it obvious that 
minimization occurs when the internal field is as much 
as possible of opposite direction and of equal magnitude 
compared to the external field; hence diamagnetism. To 
calculate it we note that doing the square gives three 
terms. The first of these give (JUS]) above. To calculate 
the other two terms we assume that the magnetic field 
is constant B c = B c e z in the region (r < R) where 
our system of charged particles resides. In order to get 
finite results we take the external field to be a dipole field 
outside some large enough radius R, so that for r > R 
we have B c (r) = (B e /2)[3(e z • e r )e r - e z ]{R/r) 3 . This 
gives us, 
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Bi ■ B c dV = —Vj ■ A e (rj), (10) 



where, 



6{r) 

A e(r) = -^-B c e z x r, 



(11) 



with 9(r) = minfl, R 3 /r 3 ), is the vector potential of our 
external field [4l|. Finally we have the constant energy 
of the external field, which is given by, 



(12) 



Assuming a constant external field extending to infinity, 
as is often done in simplified treatments, makes not only 
E e , but also -Eic, infinite. 

Summarizing, we have obtained the energy, 



E(r k , v k ) = ^ - [rrijV? + e^(rf,r k )] 



(13) 



-Ai(rj-r k ,v k ) + A e ( rj ) 



E, 



The BvL theorem is obtained by neglecting the internal 
magnetic field B\, i.e. the Aj-contribution here. 

The Darwin Lagrangian C(r k ,v k ) is obtained by 
changing the sign of the electrostatic term in (|T3)) . It 
should be emphasized that the Darwin Lagra ngian, w hich 
is presented in many advanced textbooks 1201-425 1 . de- 
scribes most of classical electromagnetism [26(, except 
that radiation and highly relativistic effects are neglected 
27]. The perfect diamagnetism of perfect conductors and 
superconductors is well described in the Darwin formal- 
ism 12, since these systems, due to the absence of 
dissipation, can be studied using classical electrodynam- 
ics. 



LEGENDRE TRANSFORM OF THE DARWIN 
LAGRANGIAN 

To investigate susceptibilities using statistical mechan- 
ics, however, requires the Hamiltonian. In principle the 
Hamiltonian % is trivially obtained from the Lagrangian 
C by means of the Legendre transform, 
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K(r k ,p k ) = Z^Pj- Vj - C(r k ,v k ), 



(14) 



after having solved for the velocities Vj = rj in terms of 
the momenta Vj(r k ,p k ), in the equations Pj = dC/dvj. 
When £ is the Darwin Lagrangian this calculation turns 
out to be difficult. Only a few texts discuss the Darwin 
Hamiltonian and mostly a first order correction, called 
the Breit term in relativistic quantum mechanics, is ar- 
rived at. The first to seriously consider the problem be- 
yond the first order approximation were Primakoff and 
Holstein (2^1. They pointed out that when v/c is not 
small there will be non-negligible effective many-body 
forces in the Hamiltonian formalism. Even if v/c is small, 
however, the terms in question can be large when many 
particles contribute, as is the case when macroscopic 
amounts of matter produce strong magnetism. Since 
strong magnetism presents no problems for the Darwin 
Lagrangian one concludes that the problem arises in the 
approximation to the Legendre transform. The present 
author has investigated this problem in detail and at least 
managed to improve the situation 3oT 35j|. 
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It turns out that one can find very good Hamiltoni- 
ans for two different limit situations. In the first case 
one finds the Hamiltonian as an expansion in the dimen- 
sionless parameter Nr e /R, where iV is the number of 
particles with correlated velocities in a region of size R 
and r e is the classical electron radius. For small values 
of this parameter, i.e. small density of charged particles 
the first couple of terms of the expansion should thus 
be excellent. The other limit that gives a definite re- 
sult is the continuum limit of a constant, not necessarily 
small, density of charged particles. Below we will con- 
sider the interaction of two systems of charged particles, 
one that produces a strong magnetic field, and one that 
responds to the energy of this field. Between them we 
assume that there is vacuum so the first version of the 
Hamiltonian, expanded to second order, should describe 
the interaction accurately. Note that this second order 
Hamiltonian, containing terms quadratic in the vector 
potential, is the one normally used for charged particles 
in an external field. Systems of many classical charged 
particles (plasmas) are still not well understood and it is 
not unlikely that the mathematical difficulties in obtain- 
ing their Hamiltonian somehow reflects this. 



positive background that make electrostatic effects neg- 
ligible. The vector field Aa of Eq. (fT6|) is then defined 
by, 

A A (r) = JL / i^l dV , (18) 
mc z J \r — r'\ 

In the simple symmetric cases that we will deal with here 
one finds that, 

— A(r) = u s v{r), (19) 
mc 

for points where g(r) =/= 0. Here v s is a dimensionless 
constant. This gives the simple result, 

A A (r) = u s A{r), (20) 

connecting the two vector fields. One notes that attempts 
to estimate v s in a concrete situation requires that the 
current density can really be seen as a product of a charge 
density and a velocity. 

THE REPULSIVE MANY-BODY FORCE 



EFFECTIVE ONE-PARTICLE HAMILTONIAN 

One of the most relevant results of these investigations 
is a qualitatively meaningful Darwin Hamiltonian which 
includes quadratic terms in the Darwin vector potential 
(note that this is a unique object which is not subject 
to gauge freedom). Since quantum investigation of dia- 
magnetism requires the square of the vector potential it 
seems natural that the traditional ("first order") Dar- 
win Hamiltonian fails to describe it. Using this improved 
Darwin Hamiltonian [30| one can determine the effective 
one particle Hamiltonian, namely all the terms that refer 
to, say, particle 1, in the the form [311 ] . 
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p--A(r) + e<t>{r) + V D (r) (15) 
c 



where cf> is the electrostatic potential energy and, 



Vb(r) - — 
mc 

Here, for our purposes, 
A(r) = 



A A 



1 



-A 



J \r-r>r ' 



(16) 



(17) 



where the current density j(r) = g(r)v(r) of all the other 
particles is the product of the charge density g with its 
velocity v. From now on we assume that we are deal- 
ing with electrons since all effects we are discussing are 
much smaller for nuclei and positive ions which conse- 
quently are assumed to simply provide a neutralizing 



To find a value for v a we must consider a specific situ- 
ation. Assume that we have a rotating superconducting 
sphere. On such a sphere surface current is induced that 
produces a magnetic field: the London moment. It is 
constant inside the sphere and a dipole outside so that 
the vector potential is given by HI]). Assuming that the 
velocity of the charge carriers is due to the rotation one 
can estimate the relevant surface charge density [28| and 
it turns out that v s is order of magnitude 1 . For simplic- 
ity we use this value below, though strictly speaking it 
should be a bit smaller for convergence of the Hamilto- 
nian inside the source. In most situations it is unfortu- 
nately not obvious how to find the amount of charge and 
corresponding velocity that produces a given j, or B. 
Nor is it clear how this is done when the magnetic field is 
due to ordered spins. In any case there is some reason to 
believe that v s is of order of magnitude unity in macro- 
scopic systems. Qualitatively a small charge density and 
large velocity makes v s smaller for a given j. 

If we use cylindrical coordinates p,tp,z, with r 2 = p 2 + 
z 2 , and assume v s — 1 the diamagnetic potential energy 
Vd of (fT6]) is, according to (fTTj). 



T/ / n _ e2 3 R 2 2 / 1 for < r < R 
Vd(V) " ^8 9 X \ i? 6 /r 6 for R < r < oo 



(21) 



i.e. harmonic attraction towards the z-axis for r < R 
and a repulsion outside the sphere, except on the z-axis 
where there is no force. 

Here we must carefully point out that the vector po- 
tential of the Darwin formalism has no gauge freedom. It 
must be divergence free, respect the symmetries of the 
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system, have a singularity at a point charge, and go to 
zero far from its sources. We note that the extra (dia- 
magnetic) term in the Hamiltonian (|15[) is proportional 
to the square of the classical electron radius r = e 2 /mc 2 
since it appears also in A a- This type of force cannot be 
found simply from Maxwell's equations and the Lorentz 
force law, but would normally require considerations of 
radiation damping and the Thomson cross section. Here 
it arises from many body effects, in spite of neglect of 
radiation. 

An even more interesting specific case is the field from 
an ideal infinite solenoid. Assuming that the radius of 
the cylindrical solenoid is p — R we have that, 



Mr) 



BR 2 



for R < p < oo, 



(22) 



where B is the constant magnetic field inside the solenoid. 
As is well known the field outside is zero. With 
our Hamiltonian (|15j) we thus find a repulsive force 
—dVojdp ~ l/p 3 , where Vd is given by, 



Vb(r) 



2 mc 



BR' 
~2p~ 



(23) 



outside the solenoid (R < p). This classical diamagnctic 
effect ~ A 2 in a field free region is different from the 
quantum mechanical Aharonov-Bohm effect, which pre- 
dicts a phase shift linear in A. Boyer [36| previously in- 
vestigated the classical problem but did not find the force 
predicted here. It is quite difficult to correctly include 
all contributions to the energy unless one starts from the 
fundamental Lagrangian and take care not to lose them. 
Nevertheless, an experimental test seems much needed. 

Let us investigate the possibility to find the force F p = 
—dVo/dp, with Vd given by (J23J) experimentally. We 
find that. 



emanating from magnetized plasmas. Recent work 37 1 
indicates that the parameter Nr c /R, limiting the valid- 
ity of our Hamiltonian, is related to the maximum elec- 
tron density (Greenwald density) above which laboratory 
plasmas tend to disrupt. At these densities the mag- 
netic energy per electron can become comparable to the 
rest energy mc 2 of the electrons. The fact that plasmas 
in general tend to be diamagnetic was noted already by 
Alfven [HI and lends support to the present results. 



CONCLUSION 

The diamagnetism found here for a system in a con- 
stant magnetic field is very similar to the quantum me- 
chanical one for a zero angular momentum atom. Since 
that quantum mechanical effect is closely related to the 
classical Larmor's theorem predicting a rotation of a sys- 
tem in an external field |16{ it is not really surprising 
that it also shows up in a classical phase space formalism 
that is careful to introduce canonical momenta systemat- 
ically. This is after all a first step to quantum mechanics. 
Diamagnetism can be powerful enough to levitate macro- 
scopic objects [39( and this indicates that some form of 
classical understanding should be possible. Pending ex- 
perimental verification of the force — VVd, and a full un- 
derstanding of the constant v s , our findings remain pro- 
visional and further research is desirable. But the main 
message here, that a Hamiltonian for a classical system of 
point charges that takes the Darwin velocity dependent 
magnetic interactions into account does predict classical 
diamagnetism, is hopefully now beyond dispute. 
Note added in proof: A possible experimental verification 
of a classical diamagnetic current can be found in 40] . 



F p (p) = 6irr e R 
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where r c = 2.82 ■ 10~ 15 m, and the term B 2 /8n in the 
first parenthesis is the magnetic energy density. In terms 
of Si-units this energy density is written, 
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(25) 



where N t is the number of turns of current / in the 
solenoid of length L. The Si-unit of this is J/m 3 . Since 
the classical electron radius and the radius R of the 
solenoid both are lengths with Si-unit meter (m) we see 
that the unit of the force is indeed J/m = N. With 
R = 0.1m, B = 15 T one finds the force F p (0.1 m) = 
4.76 • 10 _7 N on one electron. For just 10 6 classical elec- 
trons jUJ the force should thus be roughly half a Newton. 

It is tempting to see the repulsive force discussed here 
as a possible explanation of the ubiquitous stellar winds 
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